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Set Books 

D. L. Kreider, R. G. Kuller, D. R. Ostberg and F. W. Perkins, An Intro- 
duction to Linear Analysis (Addison-Wesley, 1966). 

E. D. Nering, Linear Algebra and Matrix Theory (John Wiley, 1970). 

It is essential to have these books; the course is based on them and will 
not make sense without them. 

Conventions 


Before working through. this correspondence text make sure you have 
read A Guide to the Linear Mathematics Course. Of the typographical 
conventions given in the Guide the following are the most important. 


The set books are referred to as: 


K for An Introduction to Linear Analysis 
N for Linear Algebra and Matrix Theory 


All starred items in the summaries are examinable. 


References to the Open University Mathematics Foundation Course 
Units (The Open University Press, 1971) take the form Unit M100 3, 
Operations and Morphisms. 


32.0 INTRODUCTION 


In Unit 23, The Wave Equation, we have seen how to use the one-dimensional 
wave equation in modelling the vibrations of mechanical systems such as 
the air in a musical instrument. The wave equation is only one of a variety 
of partial differential equations which are of great value in mathematical 
modelling. The intention of this unit is to introduce a few of these other 
equations. As a case study for the use of these equations in modelling, we 
shall take the conduction of heat in a solid. We shall see how different 
physical situations lead to different equations (although these are all 
special cases of one very general equation), and how equations that are 
superficially very similar may have very different properties. 


We shall begin by showing how physical problems involving heat conduc- 
tion can be modelled by partial differential equations with suitable 
boundary conditions. We shall then show how some of these boundary- 
value problems can be solved using extensions of the technique of 
separation of variables treated in Unit 23. Finally we shall compare the 
various equations with each other and with the one-dimensional wave 
equation, and indicate some of the ways in which the theory generalizes 
to more complicated cases. 


We shall concentrate on three equations: 
the one-dimensional heat conduction equation 


Ou _ 2 Ou 

at” a 
the two-dimensional Laplace equation 

Ou $ au 

ax?" dy? | — 
and the version of this equation used when the rectangular coordinates 
x, y in the plane are replaced by polar coordinates r, 0 (see sub-section 
21.1.2 of Unit M100 27, Complex Numbers I) 

au Lau bus 

or? rôðr raoe 


0 


0. 


32.1 THE FORMULATION OF PARTIAL 
DIFFERENTIAL EQUATIONS 


32.1.1 The One-dimensional Heat Conduction Equation 


Two of the principal steps in setting up a mathematical model are: 
(i) specify the unknown quantities; 
Gi) formulate the scientific laws connecting them. 


In problems concerning the conduction of heat through a solid the 
unknown quantity we wish to calculate is the temperature. It will depend 
in general on position and time and will therefore be a function u of three 
space variables x, y, z and one time variable f 


u: (x, y, z, f) ——95 u(x, y, 2, t). 


In this unit, however, we shall simplify the discussion by confining it to 
cases where the temperature depends only on two of these variables, 
either x and t or x and y, so that 


u: (x, t) —> u(x, t) 
or 
u: (x, y) —— u(x, y). 


The scientific laws relating to heat conduction are considered in the next 
reading passage, where they are used to formulate the differential equation 
of heat conduction for a one-dimensional problem. 


READ from line 9 on page K512 to line 5 on page K513; then 
READ from line 8 on page K514 to the end of the section. 


Notes 


(i) line 14, page K512 “constant on each cross section”. That is, the temperature 
in the rod is independent of y and z and can therefore be described by a function 


w: (x, t) — 5 u(x, t) 


temperature = u(x,t) 
all over this disk 


Gi) line 18, page K512 The “rate of change” is a distance tate of change. 

(iii) Equation (13-10), page K512 This is the law of heat conduction, some- 
times called Fourier's Law, and is the first of three empirical scientific laws 
used in formulating the heat conduction equation. Here AH does not mean 
“A times H” but is a single symbol. Since this symbol is used in a different 
sense in the next equation in K, it may be clearer to re-write Equation (13-10) as 


where J means the amount of heat crossing the disk (shaded in the above diagram) 
per unit time. A positive sign for J indicates that the heat is flowing to the right. 


The number k is proportional to the area of this disk, that is to the area of 
cross-section of the bar. The constant of proportionality (the value of k for a 
he clad of unit area) is called the thermal conductivity of the material of 
e bar. 

(v) Equation (13-11), page K512 This is the second of the three empirical 
laws: the amount of heat supplied to a body equals the specific heat times the 
mass times the rise in temperature. The symbol AH now has a different meaning 
from that given to it in Equation (13-10). It stands for the rate at which heat is 
accumulating in some portion of the rod whose mass is m and whose length is 
short enough to permit us to treat the temperature as uniform (i.e. independent 
of x) along it. This approximation does not lead to any error in the final result 
(as we shall see in note (vii)). 


y 

| 

| 

J —— - > 

^ x 
— J 
= mass m 

accumulating 
heat AH in 
unit time 


(v) line —7, page K512 The symbol Ax is a single symbol, like AH. 

(vi) line —1, page K512 This equation comes from the third of the three 
empirical laws, which tells us that (in the physical situation considered here) 
heat is not created or destroyed; so, the amount accumulating in any region 
equals the amount coming in minus the amount going out during the same time 
interval: 


accumulation — input — output 


Laws of this form are called conservation laws. For our problem, the conserva- 
tion law for heat applied to the portion of the bar between x and x -+ Ax shown 
in the diagram above, gives 


AH — J(x,t) — J(x + Ax, t) 


since AH is the rate at which heat is accumulating, J(x, t) is the rate at which it 
is entering the left-hand end of the portion (at time 7), and J(x + Ax, t) is the 
rate at which it is leaving the right-hand end. 


Combining the above equation with the one in note (iii), we obtain 
ðu ðu 

AH=k g (et Ax, t) - ks n) 
which is the equation at the bottom of page K512, expressed in the notation of 
Unit 23. 
(vii) lines 2, 3 and 4, page K513 In the notation used just above, this equation 
is 

1 
Ax 


ôu ou _ cp ou 
(eis 0-20) =H OO 


The reason for taking the limit of small Ax is that it eliminates the error due 
to the approximation mentioned in note (iv), where we assumed that the temper- 
ature at time f was the same throughout the portion. (A similar limit process is 
used in Unit 23 to derive the wave equation.) In this limit the left-hand side 
of the above equation approaches 


9 fou 
z C Gt j. 
by the definition of a partial derivative, and so we obtain line 4 on page K513. 


(viii) line 11, page K514 The symbol A stands for the x-coordinate of an' end 
of the rod. 


. 2u R r 
(ix) line 15, page K514 “insulated at A”. The vanishing of 5:05 t) in this 
case follows from Equation (13-10). 


(x) Fast line of the section, page K514 If one end of the rod is exposed to the 
air, it is a fair approximation that the rate of heat flow out of it is proportional 
to the difference in temperature between the rod and the air. Thus, to use boun- 
dary condition 3 we should take the zero of temperature measurement as the 
temperature of the surrounding air. 


In the derivation given, the rod is assumed to be slender and surrounded 
by insulating material. Would the same equation apply if the rod were not 
slender, e.g. if it were replaced by a slab whose two faces were each at uni- 
form temperatures (i.e. temperatures independent of y and z)? 


xt 


In fact, precisely the same equation applies. It does not even matter very 
much if the edges of the slab are not insulated, since they are too far away 
from most of the slab to have a significant effect. The need to consider a 
thin rod in the previous derivation was to ensure that the temperature dis- 
tribution is constant over any cross-section. For the slab, such a distribu- 
tion is ensured by maintaining each large face at uniform temperature. 


Exercises 


1. In the boundary condition 3 on page K514, which sign should the 
number 4 have at (i) the right-hand end of the rod, (ii) the left-hand 
end, (assuming x increases from left to right) to model the physical 
situation described in note (x) above? 


2. Exercise 2, page K514. 
Solutions 


1, If the temperature of the rod exceeds that of the surrounding 
air, i.e. u(x, t) > 0, then heat will flow out of it. At the right- 
hand end (case (a)) such a flow corresponds to J > 0 (where 


J is defined in note (iii)) and hence to x (4, t) « 0. Since u 
x 


and = have opposite signs at (A, t) we need A < 0 in this case. 
At the left-hand end (case (b)) heat flowing out of the rod 
is travelling to the left, so that J < 0, =~ (4, t)>0, and we 
need A > 0. 


2. This time the principle of conservation of heat must be 
modified to 


accumulation = input — output + production 
AH = kJ(x) — kJ(x + Ax) + g(x, t) Ax 


Using the first two empirical laws, we obtain 


Ou Qu ou 
Ax t — 
cp Ax a (x, t) {E (x + Ax, t) ax (x, 0) 


+ g(x, t) Ax. 
Dividing by Ax and taking the limit of small Ax then gives 


ou Ou 
cp a DES kao t) + a(x, t). 


A further division by k gives the equation in the question, 


T"! cp 
th — —- y 
with z for b, and k for a 


32.1.2 Laplace’s Equation 


To show how the discussion given in the preceding sub-section can be 
generalized, we now derive an equation describing heat flow in two 
dimensions. A physical situation leading to such a problem is a metal 
plate of any shape, say a circle or rectangle, sandwiched between two 
layers of insulating material, so that heat does not flow perpendicular 
to the faces, and the temperature is therefore independent of the co- 
ordinate in this direction, but depends on the other two. In the situation 
depicted in the diagram the temperature depends on x and y, but not on z. 


+— insulator 
*—- metal 


4— insulator 


To keep the number of variables down to two, we shall assume in our 
derivation that the temperature does not depend on time. (In other 
words the slab has been heated up in some way and has then been allowed 
to settle down to some “steady” temperature distribution.) The tempera- 
ture function then has the form 


u: (x, y) ——9 u(x, y) 


with domain corresponding to the region of the xy-plane occupied by the 
plate. Such a temperature function, independent of the time, is said to 
describe a steady-state situation. 


To formulate the equation we use the same empirical laws as before, but 
instead of applying them to a portion of a rod with length Ax we now take 
a rectangular portion of the plate. 


y 


xt 


(i) Fourier’s law relates the heat flows across the sides of the portion 
to the derivatives of u(x, y). Consider the side AB. An approximation to 


een ee NE i . ðu 
the temperature derivative in the x-direction at points along AB is x 


evaluated at the mid-point of AB. If we write J,, for the heat crossing 
side AB to the right, we then have 
Qu 
x -=-k>— 
Jas gx + 2 Ay) 


where k is a constant of proportionality analogous to the one used in the 
preceding sub-section (Equation (13-10), page K512). The heat crossing 
Qu , " A 
AB for a given = is proportional to the length of AB; so k will also be 
proportional to this length, that is to Ay. Calling the constant of pro- 
portionality k,, we thus have 
Qu 
Jag = =k, dy x (x, y + 3Ay) 
In a similar way the heat flow to the right across CD is 
Qu 
Jeo S =k; Ay & + Ax, y + Ay) 
For the flows across AD and BC we take the direction of increasing y as 
positive, obtaining 


Qu 
Japz -k Ax gy & + 15m) 


Qu 
Jac & —k, Ax gy (+ thx y+ Ay) 


(ii) The specific heat law is Equation (13-11) on page K512. Since we are 
considering a steady state, this reduces to the statement that heat does not 
accumulate in the rectangular portion we are considering. 

(iii) The conservation law. Since there is no accumulation, this reduces to 


0 = input — output 
= Vas + Jan) — (Jac + Jen) 
= (Jas — Jen) + (Jan ~ Jac) 


ĝu ĝu 
xk to(— Fur + n + Boot b) 


Ou Qu 
+k, mh- etas +B ors asy a) 


We now take the limits of small Ax and Ay, first dividing both sides of the 
equation of Ax x Ay so as to get a non-trivial result. 


ĝu ðu 
met Ax, y + FAY) — = 65 y + 3A) 


0= limk, 
feng Ax 
Qu Ou 
x (x + Ax, y + Ay) — = (x + 4Ax, y) 
+ lim k|? oy 
pow] Ay 


For the first term, taking the limit for Ay first (x and y are independent 
of each other) and then the limit for Ax gives, by the definition of a partial 
derivative, 


ky x (Flee y). 


A similar argument for the second term yields 


ð (du 8 [Qu 
0 - x, (2 (£4) . 2. (* 
( (3) t (5) 


which is equivalent to 


Ou au 
at y7 0. (1) 


Thus 


This is called Laplace’s equation in two dimensions. 


When we come to solve Laplace’s equation (1) by separation of variables, 
we shall find that it can be done only if the plate is rectangular. This is 
because the method only works if the boundary conditions separate as 
well as the differential equation. For this to happen values of x on the 
boundary must be independent of y, and vice versa. The only way of 
extending the method of separating variables to other shapes is to use 
other coordinate systems (where possible). The polar coordinate system, 
which we met in sub-section 27.1.2 of Unit M100 27, Complex Numbers I, 
is useful, for example, in heat conduction problems in pipes of circular 
cross-section. The polar coordinates (r, 8) of a point in the plane are related 
to its rectangular coordinates (x, y) by 


x=rcos 0 


y-rsin8 


x x 


One way of obtaining the partial differential equation of heat conduction 
in polar coordinates is to work directly from the above relations. You will 
find the calculation on pages K546 to 547, but the manipulations are rather 
heavy and of a type that may be unfamiliar to you. An alternative deriva- 
tion is to work from first principles, as in the preceding sub-section, using 


ii 


the new coordinate system. This derivation is given in the Appendix and 

yields 

16 / Ov 1 Pu 
( s) Ba 

or 


Ou lu | #u 
eei he serre eel 0 (2) 
antra Sa i 


which is called Laplace’s equation in polar co-ordinates.* 


Exercises 


l. Generalize the discussion given above for deriving Equation (1) to the 
case where the temperature may depend on time, /, as well as x and y. 


2. Write down the boundary conditions on u(x, y) that model the 
following physical boundary conditions at the edges of a metal plate 
sandwiched between two insulators in the z-direction: 


open to air 


c 


temp=0 temp=100 


insulated 


Take the temperature of the air to be 0. If you introduce any constant 
of proportionality, state its sign. 


3. Formulate a differential equation satisfied by the pressure of oil ina 
porous rock under the following assumptions. 


(a) The problem is two-dimensional. 


(b) The quantity of oil crossing unit length per unit time is pro- 
portional to the partial derivative of the pressure function with 
respect to a coordinate measuring distance at right angles to the 
unit length. 


(c) Oilis neither created nor destroyed, i.e. 
input — output — accumulation — 0. 
Solutions 


l. The argument given in the text applies with the following 
modifications. 


(a) The temperature function is now 
u: (x, y, t) —> u(x, y, t). 


(b) Change u(x, y + 1Ay) to u(x, y+ 4Ay, t), etc. in the 
equations for J4,, etc. 


* If you found the derivation of Laplace's Equation in Cartesian coordinates difficult, you 
may find the Appendix very helpful since the same method is employed. i 
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(c) The specific heat law gives 
ĝi 
AH = cp, Ax Ay e 
ôt 
where AH is the heat accumulating per unit area per 
unit time and p, is the mass of the plate per unit area. 
(d) In the conservation law, replace the zero in the left-hand 
i ĝi 
side by cp, Ax Ay x The final result is then 
au " Pu E ðu 
ax? 0y t 


where 


This is the two-dimensional heat conduction equation. 
Notice that for steady state problems it reduces to 


n : i Ou 
Laplace’s Equation, since in those cases T 0, and 


2u 

=z =0, it reduces 
ay 

to the one-dimensional heat-conduction equation. 


that if u is independent of y, i.e. 


2. For the edge on the y-axis: u(0, y) 2 0 
For the edge on x =a: u(a, y) = 100 Oel, b) 


For the edge on x-axis: oe (x, 0) 20 

ay 

J (xe [0, a]) 
For the edge on y = b: " (x, b) = hus, B) 


er Qu 
where h is a negative constant. In the third line we use y 


ð 
not z as on page K514, because the heat flow that the 
X 
insulation fixes at zero is in the y-direction. The fourth line 
corresponds to condition 3 on page K514, and the sign of ^ 
is as discussed in Exercise 1 of sub-section 32.1.1. 


3. The figure shows a rectangular portion ofthe rock with sides 
of length Ax, Ay, where x and y are Cartesian coordinates. 


Ay 
y+Ay eux B C 
y Al ID 
| | 
\ | 
| | " 
x X+AX x 
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Assumption (b) applied to the rectangle shown tells us that 
the flow 7,5 to the right across side AB is given by 


Gj 
Jag —k Ay Z (x, y + 4Ay) 
ôx 
where p is the pressure function. Similarly 


p 
Jap ~k Axa (+ dx y) 
ap 
Jez =k Ax gy Qe Mos y + Ay) 


ój 
Jen = —k Ay 2 (x + Ax, y My). 


Assumption (c) gives 
0 — 45 + Jan — Jec — Jen- 


By the same type of manipulation as in the last few lines of 
the text (but with p replacing wu) we deduce that p satisfies 
Laplace's equation. 


32.1.3 Summary of Section 32.1 


The main new idea of mathematical modelling introduced here is that of a 
conservation law. This is a law stating that same physical quantity is 
neither created nor destroyed and can be put in mathematical form by 
selecting a well-defined region of space and applying to itthe rule 


accumulation = input-output (1) 


In formulating the heat conduction equation in one dimension we apply 
this to the section between x and x +Ax; in two-dimensions with rect- 
angular coordinates, to a rectangular section of sides Ax and Ay, and so 
on. We also use Fourier’s law (heat flow is proportional to space rate of 
change of temperature) and the law that the rate of temperature rise of 
some portion of matter is proportional to the rate at which heat accumu- 
lates in it. The differential equation is obtained by dividing (1) by Ax (or 
AxAy) combining with the other two laws, and taking the limit of small 
Ax (or Ax and Ay). 


In this section we defined the terms 


conservation laws (page C7) 
steady-state situation (page C9) 
Technique 


Derive the partial differential equation describing the flow of heat in a 


body by applying Fourier’s law, the specific heat law and the appropriate 
conservation law. ` 
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32.2 SOLVING THE EQUATIONS 
32.2.0 Introduction 


In this section we consider some methods for solving boundary-value 
problems for partial differential equations of the types derived in sub- 
section 32.1.2. These methods are mainly based on the technique of 
separation of variables, which we studied in Unit 23. 


322.1 One-dimensional Heat Conduction: Homogeneous 
End-point Conditions 


The term "homogeneous boundary conditions" as applied to ordinary 
differential equations was introduced in Unit 25, Boundary-value Problems ` 
(page K458). In the context of heat conduction, homogeneous boundary 
conditions are those where the temperature at some spatial boundary, 
or the heat flow across it, or a linear combination of the two, is required 
to be equal to zero. On the other hand, if the temperature is required to 
have some value other than zero, we have a nonhomogeneous boundary 
condition. The homogeneous case can be dealt with by a straight-forward 
application of the method of separation of variables. 


One such problem is that which models a rod whose ends are both main- 
tained at temperature 0 and whose temperature distribution at time 0 
Is a given function f: 


solve 

8? 8i 

Rate x (Delo L] x D, o) 
with 

u(0,1)=u(L,t)=0 — (te(0, o) 
and 


u(x,0)=f(x) — (xel[0 Lp. 


We have already seen how to solve this problem by separation of variables 
in sub-section 23.3.3 on Unit 23. You may like to revise this solution, or 
the very similar one for the wave equation given in Section 13-4 of K, 
before proceeding to the reading passage where a slightly more compli- 
cated problem of the same type is considered. 


READ Section 13-6 from its start on page K528 to “... solution of the 
given problem” in line — 10, page K530. 


Notes 


(i) Equations (13-38), middle line, page K529 Comparing this with line 14 
on page K514 and our remarks on Exercise 1 of sub-section 32.1.1, we see that 
if the end of the rod is open to the air, h is positive (and zero if it is insulated). 
Gi) line 13, page K530 This formula for A, depends on the fact that the 
eigenfunctions of a Sturm-Liouville system are orthogonal in the appropriate 
Euclidean space (here C[0, LJ). See the equation two lines above Equations 
(12-31) on page K484. 

Gii) Equations (13-42), page K530 In writing this we take advantage of the 
fact that the eigenfunctions of the Sturm-Liouville problem form a basis (since 
the boundary conditions are unmixed). See Theorem 12-7 on page K487. 


Exercises 

1. Exercise 2, page K531. 

2. Exercise 5, page K531. The end-point conditions should be read as 
u0, 2) — uL, t) 0 — (te(0, co), 
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excluding the instant ¢ = 0, since the given boundary condition 


u(x, 


tells us that 


k 
u,(x, 0) = + cos I 


and so 


.ORX 
0) = ksin7- 


TX 


u,(0, 0) = u,(Z, 0) = 0 


is untrue. 


Solutions 


1. 
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Separation of variables with u(x, t) = X(x)T(t) gives the pair 
of differential equations (13-40) on page K529 with 


X(0) = X(Z) =0 
whence 
v 1 
XQ) = sin, a=-" @@=1,2,...) 
nrt 
TO= 4, exp (- 775) 


where n is any non-zero integer and A, is any real number. 
Since 
T(0) = A, 


the initial condition 
nux 
, 0) = k sin — 
u(x, 0) sin L 


can be satisfied by taking 1 = 1 and A, = k, and we obtain 
2 


u(x, t) =k sin G) exp (- zo) 
The separation of variables method gives the pair of differ- 
ential equations (13-40) with 

X'(0) = X'(L) = 0. 
The eigenfunctions ¥ are 


nx ant 
X(x) = cos 77, a= @=0, 1,2,...) 


nrt 
T(t) = A, exp (- Da) 
The series solution is therefore 


2 nax nt 
u(x, t) 2 }, A, — ———.]. 
(x, t) 2, cos ^y exp ( I) 


To satisfy the initial conditions we want 


co 
ksin = = Y. A, cos ZZ (x € [0, L)). 
n=0 


L 


By the theory in Sections 9-5 and 9-6 of K, which we studied 
in Unit 22, Fourier Series, the coefficients are given by 


L 
^ij k sin TE cos“ dx 
if n=1,2,... and a similar formula without the factor 2 
if n = 0. Carrying out the integrations and substituting in the 
above formula for u(x,t), we obtain the answer given on 
page K758 (in which 2m is written in place of n, since only 
even values of contribute to the sum). 


32.2.2 One-dimensional Heat Conduction: Nonhomogeneous 
End-point Conditions : 


A typical pair of nonhomogeneous end-point conditions for the one- 
dimensional heat conduction equation 


Pu 1 Qu 
SE ae (DELL) x (0, o) 
is 
0, = 
pozo CED o) a) 


describing a rod (or slab) where two ends (or faces) are maintained at 
different temperatures. 


The difference between this type of end-point condition and the homo- 
geneous type can be seen by looking for the steady-state solution. If 
u(x, t) is independent of time, it has the form 


u(x, t) = X(x) 


where X is some function from [0, L] to R, and the heat conduction equa- 
tion for u then implies 


X"(x) 2 0. 
Integrating twice, we obtain first X'(x) — c, and then 
X(x) = ex + c; 


where c, and c; are arbitrary constants of integration. Choosing these 
constants to fit the end-point conditions (1) we obtain as the steady-state 
solution 


x 
u(x,t) = 1o(1 - Üj 
If we had done a similar calculation for the homogeneous boundary 


conditions 


u(0, t) 20 
te pao} Cel Q 


we would have obtained 
u(x, t) =0 


as our steady-state solution. Thus a major distinction between homogene- 
ous and nonhomogeneous boundary conditions is that for the latter the 
steady-state solution cannot be 0. 


The differences between the two cases can also be stated in algebraic 
terms. For the homogeneous end-point conditions (2) the set of all solu- 
tions u of the heat conduction equation (on the given domain [0, £] x 
[0, co)) which also satisfy these conditions is a vector space, and we can 
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LM 32.2.1/32.2.2 


therefore proceed by looking for a basis of this vector space; the solutions 
obtained by separating variables form such a basis. For the nonhomoge- 
neous end-point conditions (1), on the other hand, the set of solutions of 
the heat conduction equation satisfying them is not a vector space, but a 
linear manifold (see sub-section 15.1.2 of Unit 15, Affine Geometry and 
Convex Cones.) The verification of this fact is the subject of Exercise 1 
of this sub-section. To specify this linear manifold we proceed in the usual 
way, expressing wu as the sum of a particular solution and a general element 
of the kernel of the associated homogeneous problem. We already have 
one particular solution, the steady-state solution, and so we can obtain 
the general solution by writing 


u(x,t) = 100( z B +(x, 1). 


It is then easily verified that v satisfies the one-dimensional heat equation 


au Ou 
a 
and the homogeneous end-point conditions 
0, ) 20 
AE peo} (elo co) 


so that the possible functions v do form a vector space and can be found by 
the method described in the preceding sub-section. 


The next reading passage gives an example of this procedure. 


READ the rest of Section 13-6, from line —9, page K530. 


Notes 


(i) line —2, page K530 It is not really necessary to use separation of variables 
for this particular solution; a more direct procedure is to set out directly to 
find the steady-state solution by writing u(x, t) = X(x). 

Gi) line 1, page K531 There is no choice about the value of À; only by taking 
À — 0 is it possible to have T a constant and so satisfy the nonhomogeneous 
boundary condition at x — 0 for all /. As indicated in the preceding note, the 
resulting solution wo is the steady-state solution. 


Exercises 


l. Verify that the solutions of the one-dimensional heat conduction 
equation, each of which also satisfies the boundary conditions 


TUI") eme 


for some fixed L form a linear manifold. 


2. A slab of thickness L (a section of the wall of a kettle, for example) 
is brought to a uniform temperature 100, and then at time 0 a block 
of ice is pushed against the face x=0, reducing the temperature 
there to 0 for all times ¢ > 0. The face x = L is maintained at tempera- 
ture 100 by boiling water. Find the temperature distribution in the 
slab at time /. Use the result: 

Í LO nnX qt 
x sin —— dx = — — cos mx. 
o L nn 
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Solutions 


l. Letu, uz be solutions of the one-dimensional heat conduction 
equation both of which satisfy the conditions 


u(0, t)= 0 


u(L,t)= 100} CED 0» 


Then if @,, &; are real numbers such that o + æ; = 1, the 
function u, defined by 
Ug = Oyu, + 05 U2 


is a solution of the heat conduction equation (since «3 is a 
linear combination of u, and uz) and it also satisfies 


u,(0, t)=0 
u3(L, t) = a4 x 100 + œz x 100 = 100 


Thus the solution set specified in the question is closed under 
the taking of affine combinations and is therefore a linear 
manifold. 


2. The boundary conditions are 


u0, = 0 
u(L, Pm (t € (0, 0) 


u(x, 0) = 100 (x e [0, L}) 
As explained in this sub-section, we start by writing 


u(x, 1) = 100(1 = i + (x, f) 


where v satisfies the one-dimensional heat conduction equa- 
tion and the boundary conditions 


(0, t) 20 
dn Zo) (te (0, o) 
v(x, 0) = 100 I (x e (0, D) 


The separation of variables calculation is the same as in 
Exercise 1 of sub-section 32.2.1 and gives 


o . [nnx n^z?t 
vs, D) = 5, A, sin 7) -Fr 


The initial conditions require 


10ž = ŽA sin (=) (x e (0, L) 
and so 
A= E [io jim (=) dx 
AE Zen _ —200(— 1) 
LJ\L n nu 
giving 


u(x, t) = 109(1 = z) 


_yyri -nnt 
asl 1) sin (2) exp ( a ) 


mam R 


(This solution has the disconcerting feature that it makes 
u(0, 0) and u(L, 0) differ by 100, whereas the problem as set 
Tequires them to be equal. The difficulty arises because the 
initial temperature distribution u(x, 0) = 100 does not satisfy 
the given endpoint conditions. This does not affect the validity 
of the solution for any other pairs (x, £t) except (0, 0) and 
(Ł, 0).) 


"Further practice: Exercises 9, 11 on page K532. 


32.2.3 Laplace's Equation: Rectangular Region 
Boundary-value problems for Laplace's equation 


Ou us 0 

ax? ss ay? 
are somewhat more complicated than those for the one-dimensional 
heat conduction equation. This is partly because the equation involves 
second-order partial derivatives with respect to both its variables x and y 
instead of only one, and partly because the boundary conditions, or even 
the shape of the boundary in the xy-plane, may be complicated. In this 
sub-section we consider a type of boundary-value problem for Laplace’s 
equation where these difficulties can be overcome by extending the methods 
we have been using for one-dimensional heat conduction problems. 


READ Section 14-2, which starts on page K548. 


Notes 


() Equations 14-9, page K548 This problem is illustrated by Figure 13-14 
on page K532. Without the condition m(L, y) — 0, the boundary conditions 
would be exactly the same as for heat conduction in a bar with ends at x — 0 
and x — L, both held at zero temperature—a problem with homogeneous end- 
point conditions, which we saw how to solve in sub-section 32.2.1. The method 
of solution here is very nearly identical. 

(i) line —14, page K550 “Fourier series... [—-L,L]" means “Fourier 
Sine series expansion on [0, L]”. 

(ii) Last paragraph of the section, pages K550-1 We find üs by the method 
described in the reading passage, with fs written for f. We find ws by solving 
Laplace's equation with the boundary conditions 


u(0, y) 0 ` u(L, y) - 0 
u(x, M) = fix) u(x, 0) =0 


by a method closely following that in the reading passage until just after Equa- 
tion (14-12) where we take Y(0) — 0 in place of Y(M) — 0. The effect of this is 


^ n à 
to replace sinh (z (M— ») by sinh (2) throughout the rest of the calcu- 


lation. We find x, and u by the same method as x, and us but with x and y, 
L and M, X and Y interchanged. 
Exercises 

l. Exercise 1, page K551. 


Use the result: 


3 


L 
L ——— 
[ 2x(x — L) sin IE dx = nn? 
0 L 


if n odd 


0 if n even 


2. Exercise 5, page K551. 
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Solutions 


ls 


2, 


As explained at the end of the reading passage, the solution 


Is U3 + u4 where u, and u4 are solutions of Laplace’s equation 
satisfying 


(a) u0, y) = ug (ZL, y) = u(x, M)=0 u(x, 0) = 2x(x — L) 
(O) u4(0, y) = ua(L, y) = ua(x, 0) = 0 u(x, M) = 2x(x — L) 


Case (a) is just that discussed in K and its solution is obtained 
by substituting the integral given in the question into Equation 
(14-15) on page K550. We obtain 


u(x, y) 
—16L? sin (nmx/L) sinh [wn(M — y)/L] 
TU quj sigo. n? sinh (nx MIL) 


For case (b) we proceed as K does until Equation (14-12) 
on page K549. We now argue as follows: since Y(0) — 0, 
we have C, — 0 (recall that cosh 0 = 1, sinh 0 = 0) so that 


Y.(y) = B, sinh (nny/L). 
To obtain the general solution we form the series 


2 . mmx . | nny 
ug(x, y) = 24s B, sin T sinh T 
We may write a, for the product 4, B, and so the boundary 
condition on y = M gives 


E 
M 
us(x, M) = J, a, sin DIY sinh ZE = 2x(x — L). 
PES L L 
Using again the integral given in the question, we obtain 
2 [-* . AMX 
a, sinh TM -I [ 2x(x — L) sin = dx 
—16L? , 
OS if n odd 


0 if n even, 
and so the series for wu, is 


ua(x, y) 


16L? sin (nnx/L) sinh (nzy/L) 
m? y=1,5/5,.... m sinh (an M/L) — 


The complete solution is 


u(x, y) = un, y) + Gs y) 


16L? ( 7 (=) 
= En pe 
L2 m» L 


sinh (nmy/L) + sinh [nx(M — m : 
i n? sinh (nn M/L) 


The solution is u + tt; + us + u4, Where apart from a factor 2 
u, and u4 have already been calculated in solving Exercise 1, 
and 14, uz are the solutions of Laplace’s equation that satisfy 


u,(0, y) = sin z, u(x, 0) = m(x, M) = u (L, 3) =0 


and 


ux(L, y) = sin Y u(x, 0) = u(x, M) = u,(0, y) =0. 


For u, the method of separation of variables gives u, = XY 
with 

X"4+4X=0, Y"-AY=0 
and 

Y(0) = Y(M) =0. 
The eigenfunctions of this two-point boundary value problem 
for Y are given by 


amy, 
Y(y) =A, — th: Wel2.3, occ 
(y) , Sin w TY 


The eigenvalues are 


and the corresponding solutions of the equation for Y are 
therefore 

QU nax nux 
X(x) = B, sinh M + C, cosh gi 


The boundary conditions on u, imply X(L) — 0 and so we 
obtain (as in the derivation of Equation (14-13) on page 
K549) 


X(x) = sinh * (L-x) 


and hence 
sul alit eae 
X(x) Y(y) = A, sinh M (L — x)sin M: 


A more general solution can be obtained by summing over n, 
but for our purposes this is not necessary since the fourth 
boundary condition, 
NL» 
u,(0, y) = sin w 


can only be satisfied by taking n = 1, to give 


; n . my 
sinh (s (L 3) sin M 


u(x, y) = em aL 
M 
A similar calculation for u, gives 
sinh E sin Z 
u(x, y) = x 
inh ZE 
sinh = 


and so the solution to the original problem is 


u(x, y) = u, + uy us +g 


. AY, nx P T 
=sin M (sinh M + sinh M (L- 3) 


LENT 


T n=1,3,5,.. 


x sinh (nzy/L) + sinh Im(M — pL] 
7? sinh (nx M]L) ) i 


32.2.4 Laplace’s Equation: Circular Region 


In sub-section 32.1.2 we stated Laplace’s equation in polar coordinates 
(it is derived in the Appendix), and commented on its usefulness for 
boundary-value problems involving the solution of Laplace’s equation 
in a region with a circular boundary. The equation is 


Ou lu 1 u : 
rat age C70 o 


Before considering its solution by separation of variables (this is done in 
the reading passage), let us consider what happens when the physical region 
on which u(r, 0) is to be defined includes the origin, so that the domain 
of the function w includes r = 0. Mathematically, this point is exceptional 
since the differential equation (1) is not defined for r — 0. Physically, how- 
ever, there is nothing special about this point and so we require u(r, 0) to 
behave in the same way at the origin as it does elsewhere in the physical 
region. In particular we require u to be continuous at r = 0. This condition 
restricts the solutions of the equation when the physical region includes 
the origin (compare the remark at the end of Solution 2, sub-section 25.2.1 
of Unit 25). You will see an example of this in the reading passage, where 
the function u defined by 


u(r,0)=cInkr with c5 0, 


whose domain does not include r — 0, is rejected as a solution of the 
equation. 


READ Section 14-3 from page K553 to line 8 on page K555. 


Notes 


(i) Equation (14-17), page K553 An alternative method of solution is the one 
used in the method of reduction of order, treated in Section 4-6 of K and in 

di 
sub-section 11.2.5 of Unit 11, Differential Equations III. We write z — E: and 
treat the equation as a first-order equation in z: 


dz z 
etre 
so that 
c 
z=- 
r 
i.e. 
du c 
dr r 


where c is a constant of integration. This is the same as the result of the first 
integration of the equation after Equation (14-17). The second integration 
gives u = c In r + ci, since r > 0, which is equivalent to Equation (14-18). 

(i) Equation (14-18), page K553 The solution (14-18) is only useful in cases 
where the domain of u does not include the origin, e.g. steady state heat flow 
in the annular region shown in Figure 14-3 on page K557 which models heat 
conduction between the inside and outside of a pipe. 

Gi) line —1, page K553 ‘‘single-valued”” means that the mapping 


(position in plane) ———> temperature 


is one-one or many-one, i.e. it is a function. It follows that 4 must map (r, 8) 
and (r, 0 + 27) to the same number, since they both correspond to the same 
point in the plane. This gives the periodicity condition 


u(r, 8) = u(r, 8 + 27) 


which applies to any solution referring to a circular or annular region, not merely 
to solutions that are independent of r. 
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(iv) line 16, page K554 “periodic boundary conditions" are defined as 
case 3 on page K478; in the present problem they would be 6(0) = Olr) and 
8'(0) = 8'(27). An equivalent condition is to give O the domain R instead of 
(0, 27] and require it to be periodic with period 27. The reason for this periodicity 
is explained in the preceding note. 

(v) line —8, page K554 We have not studied the Euler equation (it is treated 
in Section 4-8 of K) but all that is needed here is to observe that the solution space 
is two-dimensional (since the equation is of second order, and normal on (0, co)) 
and that a basis for the solution space is given by (r^, r7") if à z 0 and (I, Inr} 
if — 0. 

(vi) line —6, page K554 The functions r + r-" and r +> In r cannot 
have r — 0 in their domains and can therefore be used only if the domain of u 
does not include the origin. 


Exercise 


Exercise 3, page K557. (The same mathematics also applies to the problem 
of a long pipe carrying a liquid at one temperature and immersed in a 
liquid at another.) 


Solution 


The problem has rotational symmetry; so we expect u(r, 0) to 
be independent of 0. The only 0-independent solution of Laplace’s 
equation is the one given in Equation (14-18) on page K553. 
The boundary conditions then determine the constants of inte- 
gration c and k as follows: 


O=cinkR, 
100 = cIn kR; 
whence 
1 100 
k=— and c= =. 
R, [In (R2/A,)] 


The required temperature distribution is 


In (r/Rj) 


u(r, 0) — 100 in QR 


The points where the temperature is 50? are given by 


In (r/R,) 


502100 —— 17 
In (R2/R,) 


ie. In (R;/R)) = In (r/R,)?. 


Hence the points where the temperature is 50° form the circle 
(or cylinder) r = JRR, : 


32.25 Other Methods of Solution (Optional) 


The method of separation of variables is only one of the methods used in 
solving partial differential equations. There are Several others, including 
Laplace transforms, Fourier transforms, and various numerical methods. 
We have illustrated the use of Fourier transforms in Unit 31, Fourier 
Transforms. Yn this sub-section we indicate how Laplace transforms 
and numerical methods can be used to solve the type of problem discussed 
in this unit. 


To illustrate the use of Laplace transforms, consider again Exercise 2 
on page K531. (This is Exercise 1 of Sub-section 32.2.1.) The problem 
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is to find a function u satisfying the one-dimensional heat conduction 
equation 


and the conditions 
u0,1)=u(L,1)=0 — (te[0, «)) 


u(x, 0) = k sin = (x e [0, D). 


We can define the Laplace transform of u, a function of two variables 


by a device similar to that used in Unit 31, Fourier Transforms. Let f, 
be the function of one variable defined by 


f— u(x, t) (t e [0, o0)). 


For each x there is a function f,. The Laplace transform of f, is 


uo = [ea " 
o 
Since 
S(t) = u(x, t) 


we may regard (1) as defining the Laplace transform of u. We can write (1) 
as 
ES 
£Lf. Ks) = Liux, s) = f e "ux, t) dt. 
0 


We can then take the Laplace transforms of both sides of the heat con- 
duction equation, obtaining 


Ou mU 
£ [| =a : 
If we assume, for the left-hand side, that differentiation under the integral 


sign is justified, and apply Theorem 5-4 on page K187 to the right-hand 
side, we obtain 


z E[u], s) = a? Gf [uYo, s) —u(, 0). 


This equation can be written, using the initial condition on u, as 


a? 2 apa 7X 
(a —a s) otc s)=-a ksin 


and solved as an ordinary differential equation for each value of s, subject 
to the boundary conditions 
£[u](0, s) = £[u](L, s) = 0. 


The solution is 


? nx 


a : 
TMG 9 7 2 p y 


as you can check by substituting into the differential equation. Finally, 
we invert the Laplace transform using the table on page K229, and obtain 


. TX 
u(x, t) = k exp (—7^ t[L^a?) sin T 
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For this particular problem separation of variables yields the solution 
more quickly, but for more complicated problems the Laplace transform 
method often turns out to be a useful source of information about the 
solution. 


To illustrate the numerical solution of partial differential equations let us 
consider the solution of Laplace’s equation 


3u Ou 
æ tay 
on the triangular domain shown below subject to the conditions 
u(x, 0) = x?, 
u(0, y) = 0, 


u(x, 1—x) =x. 


y 


x 


Separation of variables will not work here because there is no suitable co- 
ordinate system (x,, x2) for which the boundaries can be described as x; = 
constant. For a numerical solution we choose a step size / and look for 
approximations to w(x, y) at the points where x and y are integer multiples 
ofh. 

The diagram shows these points for the case A = 4, together with those 
values of u (in red) that are determined directly by the boundary conditions. 


ty 
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The unknown values H4, Ug, Uc, Where u4 = x(4, 4), etc., can now be 
calculated approximately using an approximate formula for the deriva- 
tives in Laplace’s equation. This approximation is analogous to the 
approximation for the second derivative of a function in C?(R) 


re) _fx+h- amm fx —h) + Od?) 


which we used in the optional sub-section 21.4.2 of Unit 21, Numerical 
Solution of Differential Equations. It is 


eu u(x + h, y) — 2 m -h y 
I ye X) aiz Y) + u(x ny) oh?) 


eu (y= u(x, y + h) — 2u(x, y) + u(x, y — h) 
ay! 7" h? 


so that Laplace’s equation becomes (neglecting the terms of order l?) 


+ O(n?) 


u(x + h, y) + u(x — h, y) + u(x, y + h) + ux, y — h) — Au(x, y) = 0 
In other words, the value of u at each mesh point of a network like the 
one shown in the diagram is approximately the average of its values at 
the four neighbouring mesh points. 
Applied to the grid shown above, the approximation requires 

(at A) 0.50 + 0 + 0.25 + ug — 4u, = 0, 

(at B) uc +0 + u4 + 0.06 — 4uz = 0, 

(at C) — 0.75 + ug + 0.50 + 0.25 — 4uç = 0. 


By solving this system of equations (using, for example, Gauss elimination 
with interchanges) we obtain 


44 =0.23, ug 20.18, uc =0.42. 


The accuracy of an approximate solution obtained by this method can 
be improved by going to a smaller value of h, say half the previous one, 
and the error estimated by comparing solutions with different values of /i. 


32.2.6 Summary of Section 32.2 


For homogeneous boundary conditions (i.e. boundary conditions compat- 
ible with a zero solution) the methods described in Unit 23 were used. For 
the heat equation with nonhomogeneous boundary conditions we look for 
a solution in the form 
u=u,+tov 

where u, is the steady-state solution, i.e. one that does not depend on the 
time variable f, and v is a solution which satisfies a homogeneous boundary 
condition. For Laplace's equation in a rectangular region we look for a 
solution of the form 


u = AG +U ty a 
where u; satisfies the given boundary condition on the ith side of the rect- 
angle and is zero on the other three sides. For Laplace's equation in a 


region with a circular boundary we use polar coordinates and apply the 
boundary condition that v must be a periodic function of the angle 0. 


No new terms were defined in this section. 


"Techniques 


1. Use the separation of variables technique to solve the one-dimensional 
heat conduction equation with homogeneous and nonhomogeneous 
boundary conditions. 

2. Solve Laplace's equation for rectangular and circular regions by 
applying the separation of variables technique. 


323 PARTIAL DIFFERENTIAL EQUATIONS 
OF SECOND ORDER 


32.3.0 Introduction 


We have studied a number of different second-order partial differential 
equations, but the problems and results were somewhat different in differ- 
ent cases. In this final section we look at partial differential equations of 
second order in a more general way, to understand better the reasons for 
these differences. We begin by comparing the main equations we have 
been studying, and the boundary conditions we have used for each. 


32.3.1 Well-posed and Ill-posed Boundary-value Problems 


Written in a standardized notation to make them easier to compare 
and with a — I, our three main partial differential equations are: 


the wave equation 


Pu Pus 
Ox xl 

the heat conduction equation 
^u E Qu 0 
Ox Ox, 


Laplace's equation (steady heat conduction) 
du Pu 7 
Ox Gx2 C 
(We shall bring Laplace's equation in polar co-ordinates into the story 
in the next sub-section.) 


Superficially the equations look quite similar, but in fact they are very 
different; in particular the type of boundary condition used is different 
for each equation. Typical boundary conditions are illustrated below. 


xt 


u given 


1 u given u given 
Bu. gwen 
ax, 9! 
Wave Equation Heat Conduction Laplace's Equation 
Equation 
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For the wave equation we have two conditions when X3 — 0; for the 
heat equation, one; and for the Laplace equation we have one condition 
again but this time there is a second on the line x, = M. When we ap- 
proached the equations from physical problems, these differences arose 
naturally from the differences in the problems; but the object of the present 
section is to relate these differences to the natures of the equations them- 
selves and not to the physical situations that they model. 


We can understand the differences if we look at the solutions obtained by 
separating the variables, writing 


U(X, x2) = X(x) X(x). 
The equation for X, is, in each of the three cases, 
Xi + 2X, =0 
with solutions 
" nnx; 
X(x) =A sin C. where A is arbitrary, 


and eigenvalues 
A——— (021,2, ...). 


corresponding to the operator — D?. 


The equation for X, is different in each case. For the wave equation it is 


Xj +AX,=0 Q) 
For the heat conduction equation it is 

X,+4X,=0 (2) 
and for Laplace’s equation it is 

X;—2AX,20 (3) 


For Equations (1) and (3) we need two boundary conditions to determine 
a solution, and for Equation (2) we need one: but just what type of 
boundary conditions are most suitable? Remember that, because of 
the condition A — n^z?/L?, we must be prepared for very large positive 
values of 2. 


We have a clue to the answer from our study of ill-conditioning in Unit 7, 
Numerical Solution of Simultaneous Algebraic Equations, and Unit 21, 
Numerical Solution of Differential Equations. We saw there that if a re- 
currence relation or differential equation has an increasing exponential 
solution in its xernel, then any initial-value problem for it will be (abso- 
lutely) ill-conditioned; but if the solutions are bounded or decreasing as x 
increases, then the initial-value problem will be (absolutely) well-condi- 
tioned. Applying this criterion to Equations (1), (2), (3) above, we see 
that initial-value problems for Equations (1) and (2) should be well- 
conditioned, for all 2, since the kernel of the operator D? + 2 associated 
with Equation (1) is spanned by the bounded functions x ——— sin Jax, 
and x > cos ./Ax, and the kernel associated with Equation (2) is 
spanned by the decreasing function x ———» exp (— 4x). On the other 
hand, for Equation (3) the kernel is spanned by x ———» exp (fax) 
and x +——> exp (C 23), and since one of these is an increasing expo- 
nential we would expect an initial-value problem in this case to be ill- 
conditioned. Arguing back from Equations (1), (2), (3) to the partial 
differential equations they came from, we may expect, therefore, to find 
that initial-value problems are well-behaved for the wave and heat con- 
duction equations, but not for Laplace's equation. And it is indeed for 
the wave and heat conduction equation that we have been using initial- 
value problems (see the illustrations above), with two initial conditions 
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for the wave equation because its order in x; is two and one for the heat 
equation, whose order in x, is only one. 


This argument can be carried further. We saw in sub-section 7.3.2 of 
Unit 7 how to formulate well-conditioned boundary-value problems for 
second-order recurrence relations with one increasing and one decreasing 
exponential solution in the kernel. These were the two-point boundary- 
value problems. Just the same considerations apply to the second-order 
differential equation (3) 


X1 — AX, =0. 


If X,(0) and XY,(M) are given, we have a two-point boundary-value 
problem which is well-conditioned. Arguing back to the partial differ- 
ential equation, we may expect to find that a two-point boundary-value 
problem for Laplace's equation, with the value of u(x,, x2) specified 
for x, = 0 and for x; = M, to be well-behaved; and it is indeed for Lap- 
lace's equation that we have been using this type of boundary condition. 


The conclusions we have drawn from our knowledge of ill-conditioned and 
well-conditioned problems are confirmed by the general theory of partial 
differential equations. This theory includes results which do for partial 
differential equations what the existence and uniqueness theorem 
(Theorem 3-2, page K104) does for ordinary differential equations. In 
particular it tells us that the three boundary-value problems illustrated 
in the above diagram have unique solutions, provided the given boundary 
values satisfy suitable conditions (e.g. continuity) These boundary- 
value problems are said to be well-posed. On the other hand if we try to 
specify the boundary values in one of the ways that led to ill-conditioning 


in our discussion—for example by giving values of u(x,, 0) and x (xı, 0) 
2 


on the boundary x; — 0 for Laplace's equation instead of for the wave 
equation—then it is quite possible that the boundary-value problem has 
no solution at all. Such problems are said to be ill-posed. 


Exercises 


l. It is proposed to calculate u(x,, x2) for 0 < x, « 1,0 & x; <1, given 
that satisfies the equation 


Ou ðu 

axi + x; -0 
with 

u(0, x2) = u(1, x2) =0 (x; € [0, 1]) 
and 


4%1,0)=f@) — GG € [0, 1) 


where f is a given continuous function, with FO) — f(1) = 0. Using 
an argument based on separation of variables similar to that given 
in the text, would you expect this problem to be well- or ill-posed ? 


2. Consider the following boundary-value problem: 
u = ðu 
ax} ax} 
with 
u(0, x2) =u(1,x2)=0 (x, € [0,2]) 


ulxy, 0) = fle) 
Ln cul Gr e DD, 1) 


where f and g are given continuous functions, with FO) =f) = 
&(0) = g(1) =0. 
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(i) Aou doing any calculation, and relying on the results stated 
above, would you expect this problem to be well- or ill-posed? 


(i) Assuming that f &, try to carry through the solution of this 
problem by separation of variables. 


Solutions 


l. 


Separation of variables with 


u(x, x2) = X(x) X(x) 


gives 

Xi +4X%,=0 with X0) = ¥,(1) =0 
so that 

XQ) = Asin nax, (x, € (0, 1) 
and 


A=nn? (n-12,3,...) 
The equation for X; is 

X1—AX,-0 
so that 

X3(x5) = c exp (n?1?x;). 


Since this solution is an increasing exponential, the initial- 
value problem for X, is ill-conditioned, so we expect that 
problem of finding u(x,, x2) from its initial values at x, = 0 
to be ill-posed. There is, in general, no solution. 


() For the wave equation we should have the initial values 
of u and its derivative, not the initial and final values of u. 
The problem is ill-posed. 

(i) Separation of variables as in Solution 1 gives the same 
solution for X;, with 4 — n?z?, and the equation for 
Xis 

XZ + AX, =0 
giving 
X2(X2) = c, cos nux; + c; sin nx; 


which implies X;(0)- X,(2). Thus every solution 
obtained by separation of variables has 


u(x, 0) = u(x;, 2) 


for all values of x. The same will be true of any linear 
combination of such solutions. Thus a solution is 
possible only if f and g are the same function, contrary 
to the assumption that f# g. 


32.3.2 The Classification of Partial Differential Equations 


In our discussion of boundary conditions in the preceding sub-section 
we did not give separate consideration to Laplace’s equation for polar 
coordinates. Since this equation is merely the Cartesian form of Laplace’s 
equation transformed to a new coordinate system, we would expect this 
equation to have the same essential properties as the Cartesian form: 
that is, we would expect to find that a boundary-value problem with the 
value of u specified on a closed boundary is well-posed and that a problem 
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with both u and = specified on some part of the boundary, is not well- 
X 


posed. This turns out to be the case. In fact we can go further and specify 
a general class of equations for which giving the value of u(x,, x2) for 
(xı, x3) on the boundary gives a well-posed problem. These equations are 
called elliptic differential equations. 


The general linear homogeneous partial differential equation of second 
order is 


Ou au Ou ôu ôu 
= +2 +wu=0 
P Ox? +29 Ox, 0x; bad 6x3 à “ox, ôx 


where p, q, r, s, t, w are given functions with the same domain as u. An 
equation of this type is said to be elliptic if the quadratic form 
(E, M) ——9 Plai, x2)? + 2405, xi)&n 
*rGsx)r) (n) ER?) 


is positive definite* for all x,, x; in the domain of u (because the graph of 
such a quadratic form is an ellipse). For example, with Laplace’s equation 


Ou ou 0 
adag 


the quadratic form is 
(m — @ +? 
and for Laplace's equation in polar coordinates 


Ou 1 Qu 1 Ou 
L qnte. m 0 
Ox; x, ax, + xióxi 0 oe 
it is 


Qo — gale 
x, 

Thus both forms of Laplace's equation are elliptic. The significance of 
the quadratic form is this: by changing to a new coordinate system (e.g. 
from rectangular to polar coordinates) we can alter the form of the equa- 
tion, but all that happens to the quadratic form is a congruence transforma- 
tion (see sub-section 14.1.2 of Unit 14, Bilinear and Quadratic Forms) 
and so its positive definite character cannot alter. The proof of this fact 
is laborious, so we omit it. 


The simplest way of recognizing an elliptic partial differential equation is 
to consider the function pr — g? (which is the determinant of the matrix 
of the quadratic form). If pr — q? takes positive values for all (x, x2) in the 
domain of u, then the equation is elliptic. 


Similar considerations apply to the wave and heat equations. For the wave 
equation, the quadratic form is 


©, n) — @ — 9? 


with rank 2 and signature 0. Any partial differential equation whose 
quadratic form has rank 2 and signature 0 is said to be hyperbolic (because 
the graph of such a quadratic form is a hyperbola). Hyperbolic differential 
equations may be recognized by using the fact that if pr — q? takes negative 
values for all (x,, x2) in the domain of u, then the equation is hyperbolic. 


* Sec page N 168. 
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For the heat equation, the quadratic form is 


(¢, 1) — e 


with rank 1 and signature 1. Any partial differential equation whose 
quadratic form has rank | and signature 1 is called parabolic (because the 
graph of such a quadratic form is a parabola). Parabolic differential 
equations may be recognized by using the fact that if pr — q? has the value 
zero for all (x;, x;) in the domain of u, then the equation is parabolic. 


For hyperbolic and parabolic equations the problem of finding u inside a 
region when its values are given on the boundary is not well-posed. There 
is a general theory for these equations which specifies types of boundary- 
value problems for these equations which are well-posed, but it is beyond 
the scope of this course. 


Exercises 


l. Classify the following equations as elliptic, hyperbolic or parabolic. 


A Ou Ou Ou Qu 
i) >5+4—— sats = 
o oxi Ox, 0x; * oxi * Ox, +xju=0. 
T Pu 
(ii) pr Ax 0. 
wa Otu Ou eu Ou 
(ii) adt Lan pt gr e 


2. This question requires a familiarity with the chain rule for partial 
differentiation. If you are not familiar with this rule, treat the exercise 
as an example or omit it completely; you will not be examined on it. 

(i) Given that a function u 


u: (X1, x2) ——9 ux, x2) (Q5, x2) e R?) 
satisfies 
aus Ou rot ng 
Pad aA, Ad 
where p, q, r are real numbers, what equation must the function 
v, defined by 
u(x,, x2) = v(ax, + bx;, cx, + dx;) Qu, xe R) 


satisfy? Note that the two functions u, v represent the same physi- 
cal quantity described in two different coordinate systems. 


(Hint: 

Qu 

ax (X1, x3) = (aD, + cD3)v(ax, + bx;, cx, + dx3) 

X1 
where 
Ov 7 
Dv, y2) =g (D, ¥2) € R)) 
i 
(ii) Show that if the matrix P= [; ‘| is non-singular, the quadratic 


forms associated with the equations for u and v are congruent 
(i.e. their matrices U and V are related by V = P" UP), and hence 
that the equations for v and v are both of the same type (both 
elliptic, or both parabolic, or both hyperbolic). 


33 


LM 32.3.2 


Solutions 


Gi) 


(iii) 


2. (i) 


(i) 
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Elliptic: since the quadratic form 
(E, n) —— P + dn + 5? = (E + 20)? + n? 
has rank 2, signature 2. 


Hyperbolic: since the quadratic form 


(5n: > én = 4(č + n}? — 8E — n)? 


has rank 2, signature 0. 


Parabolic: since the quadratic form 
(5, m —— € + 46 + 4n? = (E + 2n)? 
has rank 1, signature 1. 


Since u(x,, x2) = v(ax, + bx; , cx, + dx;), we find 
Qu 
dx, Ci x2 = (aD, tcD3w...,...) 
X, 
(as given) 


Qu 
x, (x1, x2) = (bD, + dD3yC.., ...) 


eu " 

ad G1, x2) = (aD, + cD2)*0(...,...) 
Ou 

6x, Ox, Ga, x2) 

= (aD, + cD2)(bD, + dD3yw...,...) 


eu 
oe (Œr, x3) = (bD; + dD3wC..., ...) 
and hence 
2 Ou Pur 
(e ad *2q ix, Ox; tr ao) x) 
— (p(aD, + cD5Y* + 24(aD, + cD) 
X (bD, + dD;) + r(bD, + dD2)*)v(...,...). 
so that the differential equation for v is 
((pa? + 2gab + rb?) D? + 2(pac + q(ad + cb) 


+ rbd)D,D, + (pc? + ged + rd?) D3)v = 0, 


The matrix of the quadratic form associated with this 
differential equation is 


ve po + 2gab + rb? pac + q(ad + cb) + rb 
pac + q(ad + cb) + rbd pc? +2qcd + rd? 


and it is easily verified that 
V =PTUP 


where 


and U, the matrix of the quadratic form associated 
with the original differential equation is 


val? 4 
4 ry 
The matrices U and V therefore have the same signature, 
and the equations are of the same type. 


32.3.3 Summary of Section 32.3 


A boundary-value problem for a partial differential equation is said to be 
well-posed if it has a unique solution; il/-posed if it has no solution or more 
than one. The problem of solving Laplace’s equation inside a region given 
values for u on the boundary of the region is well-posed. The problem of 
solving the heat equation with u(0, f) and u(L, f) given for all t > 0 and 
u(x, 0) given for all xe [0, L] is well-posed, and so is that of solving the 
wave equation with u(0, t) and u(L, t) given for all t > 0, and u(x, 0) and 


a 
E (x, 0) given for all x e [0, £]. 


A partial differential equation of the form 


eu Ou Ou —Qu Qu 
Pat ox. * "aa uera ‘at wu =0 


(with u: (x1, x) — u(x, x2)) is said to be elliptic if pr > q?, parabolic 
if pr =q* and hyperbolic if pr « q?. Boundary-value problems for elliptic 
equations are analogous to those for Laplace’s equation; for parabolic 
equations, analogous to those for the heat equation; and for hyperbolic 
equations, analogous to those for the wave equation. 


In this section we defined the terms 


well-posed boundary-value problem (page C30) 

ill-posed boundary-value problem (page C30) 

elliptic differential equation (page C32) 

hyperbolic differential equation (page C32) 

parabolic differential equation (page C33) 
Techniques 


l. Determine whether a boundary value problem is well- or ill-posed. 
2. Classify a given partial differential equation as elliptic, hyperbolic 
or parabolic. 
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324 SUMMARY OF THE UNIT 


We first of all constructed the partial differential equation which would 
give the heat distribution in a body as a function of position and time. 
The one dimensional case resulted in an equation of the form 

9u —,0u 


L——— 
ax? at 
called the one-dimensional heat conduction equation. We then looked at a 
steady-state situation in two dimensions and derived the equation 
9u Pu " 
ox! oy — 
which is called Laplace's equation in two dimensions. However, we found 
that for problems defined on regions with circular boundaries we have 
to re-express Laplace's equation in terms of polar co-ordinates 


2 2 
E Tor LEN ed (r #0). 
ðr? rór r? 66? 
All these equations are soluble using the separation of variables technique 
and examples of all three equations were given to show how this is done. 
We also showed (in an optional sub-section) how Laplace transforms and 
numerical methods can also be used to solve these equations. 


In the last section we looked at general second-order partial differential 
equations in a general way to see how and why the different boundary 
conditions imposed produce different solutions. This led on to the question 
of whether a problem is ill- or well-posed which is dependent on how the 
boundary conditions are specified. 


Definitions 
conservation laws (page C7) * 
steady-state situation (page C9) *o* 
well-posed boundary-value problem (page C30) * 
ill-posed boundary-value problem (page C30) * 
elliptic differential equation (page C32) * 
hyperbolic differential equation (page C32) * 
parabolic differential equation (page C33) * 

Techniques 

l. Derive the partial differential equation describing the flow of heat in * 


a body by applying Fourier's law, the specific heat law and the appro- 
priate conservation law. 

2. Usethe separation of variables technique to solve the one-dimensional * 
heat conduction equation with homogeneous and nonhomogeneous 
boundary conditions. 


3. Solve Laplace’s equation for rectangular and circular regions by *ox 
applying the separation of variables technique. 

4. Determine whether a boundary-value problem is well- or ill-posed. * 

5. Classify a given partial differential equation as elliptic, hyperbolic * 


or parabolic. 
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32.5 SELF-ASSESSMENT 


Self-assessment Test 


This Self-assessment Test is designed to help you test your understanding 
of the unit. It can also be used, together with the summary of the unit, 
for revision. The answers to these questions will be found on the next 


non-facing page. We suggest that you complete the whole test before 
looking at the answers. 


1. 


In the derivation of the one-dimensional heat conduction equation 
in sub-section 32.1.1 we assumed that the thermal conductivity k is 
constant. What would the equation be for a rod composed of a 
nonhomogeneous material whose conductivity varied along its length? 


Is it still possible to solve the problem by separation of variables and 
an eigenfunction expansion? 


A disk of radius 5, initially at temperature 0 throughout, is brought 
into contact at time t=O with a heat source, which maintains the 
temperature of the rim at 1000. Find the temperature at points in the 
disk for time t » 0. Hint: The solution of 

rR' +R +2rR=0 
which is bounded for r < b is R(r) = Jo(Ar). Also 


S Jon) ; 
Ew =? for  re(0,5) 


where 4,,42,43,... are the solutions of Jo(2b) = 0. The functions 
Jo and J, are known as zeroth and first-order Bessel functions. 
Laplace's equation in threc dimensions can be written 
l 1 
„2 a in óD EN 
[D D) + ind D, (sin $D4) + d$ 
Find a series solution by the method of separation of variables, assum- 
ing that u is independent of 0. Hint: let cos $ = p. 


Da} u(r, b, 0) = 0. 


Exercise 2, page K557. 
Exercise 1, page K531. 
Exercise 10, page K532. 
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Solutions to Self-assessment Test 


1. The point at which the derivation differs from that on page K512 is 
line —1. Instead, we should write 
ĝu 


aa- -|-4z 


x Ax Ox 


where & is a function with domain R and codomain R*. It then follows 
that 


1 Qu 
k 
xl Ox 


xt Ax ax 


Taking the limit for small Ax, we obtain 
ô ĝu Qu 
x (9 ze0- oe» 


This is again a second-order linear partial differential equation. 
If u(x, t) = X(x)T(t), then Y satisfies 
(kX) 2AX 


which is clearly in self-adjoint form. With suitable boundary condi- 
tions this is a Sturm-Liouville problem and its solutions form a com- 
plete set of eigenfunctions. For example, if k: x—- 1 — x? and 
the rod occupies the region x € [—1, 1], Xisa Legendre polynomial. 


2. The heat conduction equation is best written in polar coordinates 
1 1 A 
Ur += Uy + Mag = au, (r #0). 


Since the initial temperature distribution is independent of 6, u will be 
independent of 6 at all #, i.e. ugg = 0. We must solve 


Upp + : u,=a°u,  ((rte€(0,b) x R*). 
Let us write u(r, t) = R(r)T(t). Then 
[ro + : ro] T(t) =@RNT(N). 


Thus 
rR'(r) + R'(r) + 2rR() =0 (1) 


and 
a 
T’+-—35T=0 
a 


The boundary condition is u(b, t) = 1000 which is nonhomogeneous. 
We first obtain a particular steady-state solution with this boundary 
condition. 


u(r, t) = 1000. 


To this we add the general solution with the homogeneous boundary 
condition u(b, t) = 0, i.e. R(b) 2 0. We must solve (1) subject to this 
condition. Now the general solution of (1) which is bounded for the 
disk (r < R) is Jo(Ar); in general Jo(4b) # 0. So we have an eigen- 
value problem—only those values 4, which satisfy J (4,b) — 0 are 
admissible. The corresponding solutions for T are 


t im g~ lha) 
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3. 


4. 


Thus the solution of the homogeneous problem is 
ur, 1) = Y A Jur) e 9 (Cr, t) E (0, b) x R*) 
n=1 


The solution of the nonhomogeneous problem is 
u = u 4 yo 


and we require u(r, 0) = 0, r e (0, b), i.e. 


M 
1000 + Y, 4,J,(4,r) = 0. 
n=1 


Now 
> E ni for re(0,b) 
therefore 
|. —2000 
^o AMbAQD) 
and 


e MAS) ndm 
u(r, t = 1000/1 L HAr) caver j. 
ad PERO»; 


If you are stuck, then, before reading on, consult sub-section 25.2.3 
of Unit 25. 


The general solution is 
eo 
Y (A,r" + Br 0*9P(cos $). 
n=0 


For the method, see page K559. 


Let 
Rir" (reR*) 
R: r — r7 (re R*). 
Then 
R: r nr"! 
Rizr n(n — Dr7?. 
and 


r2 RI + rR, — HRS r n(n — Yr" + nr” nr =0. 
Also 

R,: r — -nrt 
TO8-4 


R3: r —> n(n + yr 


and 


rR; + rR, — n? Ra: ri n(n + Dr" — a — nr" =0. 


Separation of variables with u(x, 1) = X(x)T(t) gives the pair of 
differential equations 


X"—ÀX-0 
Loa 
T'-T-0 
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The conditions X(0) = X(Z) = 0 imply 


—-nig? 
X(x) = sin a I (1—1,2,...) 
-rnt 
T(t) = A, exp ( Lat 


The series solution is therefore 
2,42 
e Qo ([nnx —n"n't 
u(x, t) = LA. sin ( T ) exp ( Lai ) 


The initial condition gives 


d (ex ee 
æ . nnx 2 
u(x,0) = Y A, sin — = | 
n=1 L Z 
L-x z&x«L 
2 
so that 
Al u2 — nn% Í . onnx 
A,== x sin — dx + (L- ») sin ax} 
L [ L ua L 
4L . nx 
— -——Ssin— 


4L 2 1 —n'n*t\ . (nux 
— $ exp (=r) sin Cz) (n odd) 
0 (n even) 


Since the boundary conditions are nonhomogeneous we first solve 
the steady-state problem by looking for a solution of the form 
u(x, t)= X(x). For this u to satisfy the heat conduction equation, 
X must satisfy 


X'-0 
and the boundary conditions give 
X'(0) 2 XL) =k. 


The differential equation for X gives X(x) = €,X + c; and the boundary 
conditions give c, =k. Any value of c; will do, so we take c; -0 
for our particular solution, obtaining 


up(x, t) = kx. 
The given problem now reduces to 

u(x, t) = kx + v(x, t) 
where v(x, t) satisfies the one-dimensional heat conduction equation 
with 

20, 0-2, f)20 


v(x, 0) = —kx. 
Separation of variables gives 
v(x, t) = X(x)T(t) 
with 
X'—AX and X(0)- X(L)-0 
T — AT. 
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so that 


v M i 
Xo) c sin T, imet dud 


25.2. 
T(t) = A, ex(- 7; Jl 


2 


The series solution is therefore 


= AMX wrt 
»1)= 2, A, sin| — == 
v(x, t) li sin( i ) exp( T7 ) 


The initial condition gives 


—kx = v(x, 0) = 4 sin 


so that 


2 nax 
A,= gd (- kx) sin 7 dx 


The solution is therefore 


o ED 2, 2r 
u(x, t) = kx + z PE D sn) exp(-" 7 ). 
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32.6 APPENDIX 
Derivation of Laplace’s Equation in Polar Co-ordinates 


As for the derivation of Laplace’s equation in rectangular coordinates 
(in sub-section 33.1.2), we consider a metal plate sandwiched between 
two layers of insulating material so that heat does not flow perpendicular 
to the plate. In addition we assume a steady-rate temperature distribution. 


| 


We apply the heat conservation law to a portion ABCD of the plate like 
the one shown in the diagram. This particular shape ABCD is chosen 
because it is convenient to the coordinate system, just as the rectangle is 
suitable for Cartesian co-ordinates. Applying Fourier’s Law we obtain 
for the flow across AB (in the direction of r increasing) 


[7 

Jap -k AB (t, 0 + 4 A0). 

But the length of arc AB is rA@, and so 
Ou 

Jas = — kar 807 (r, O 4-3 A6). 

Similarly we obtain for the flow across CD (in the direction of increasing r) 
Qu 
Je = —ky(r + Ar) AO= (r + Ar, 0 + 1 A0. 


To apply Fourier's law to the sides 4D and BC we note that in all cases 
upto now the heat flow across a line has been 


k, x (length of line) x (the rate of change of temperature with 
distance measured at right angles to the 
line). 


A small change ¢ in @ corresponds to a change rọ in AB. So the rate of 
change of temperature with distance perpendicular to AD is 


. Uulr,d+o)—u(r, 0) 1du 
1 = 
$50 rp 7 80 
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For the sides AD and BC Fourier's law therefore takes the form 


1 du 
Jap = —k, Ar ag (r+ 2 Ar, 8) 


1 ĝu 
Jae —k, Arzt t i^n 6 + A0). 


The second of our empirical laws tells us, as before, that for a steady state 
there is no accumulation of heat, and this, combined with the last (the 
conservation of heat) tells us that 


0 = input — output - 
= (Jas + Jan) — Vac + Yen) 
= —Jep + Jas — Jec + Jan 


s bol (r+ BF (r+ Ar, 0340) —r 52 (7,0 +4 80) 


1 fôu Ou 
+key Are [spe +H Ar ean - ene]. 


Dividing by Ar A0 and taking the limits of small Ar and A6, we obtain 
ð [ Ow 1 ĝu 
0-k, =r 2 Ie 
tar ( =) +h, r 00? 
which reduces to 
Ou lôu 1 du 
ôr? r Arr? a6? 


which is Laplace’s equation in polar coordinates. 


(r #0) 
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